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Recent detections of Fanai-off-Riley Class I AGNs by HESS, MAGIC, and VERITAS suggest 
that very-high-energy 7-rays (VHE, E>100 GeV) may not have a leptonic origin. We present a 
hadronic model to describe the TeV photons as the neutral pion decay resulting from p7 and pp 
interactions. For the p7 interaction, we assume that the target photons are produced by leptonic 
processes and apparent at the second spectral peak. For the pp interaction we consider as targets 
the thermal particle densities in the lobes. We show that this model can describe the TeV spectra 
of the radio galaxies NCG 1275, M87 and Cen A. 
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Fanaroff & Riley Class I (FRI) sources are radio loud active galactic nuclei (AGNs) exhibiting 
clear structure of a compact central source, twin large-scale jets and lobes. The spectra of FRI radio 
galaxies are generally well described by standard non-thermal Synchrotron and Synchrotron-Self 
Compton (SSC) models [|], The radio through optical emission originates from synchrotron 
radiation while the X-ray through y-ray emission originates from SSC. A leptonic model with only 
one population of electrons predicts a spectral energy distribution (SED) that can not extend to 
> 10^^ Hz 1^ |5|]. Also, some authors have suggested that the GeV to TeV emission may have 
origins in different physical processes [^. Therefore, we use these two facts and consider hadronic 
processes (pp and py) to describe the VHE spectra of Cen A [|, |[ |], M87 [|, ^ and NGC1275||l^, 
11]]. We first require a description of the SED up to MeV/GeV [|, |, |[ (T^. Then, we assume the 



thermal particle densities in the lobes and the lobe's distance to the core [12, 13] to describe the 
TeV energy contribution to the spectrum either by p 7 or pp interaction. 

NGC 1275 is located at the centre of the Perseus cluster at a redshift of z = 0.0179. It has a 
strong, compact nucleus and jet. Its jet has been detected with an intrinsic velocity of 0.3c — 0.5c 
and oriented at an angle of 30° - 55° to the line of sight [p^]. M87 is located in the Virgo 
cluster of galaxies at a distance of ~ 16 Mpc (z=0.0043). It has been detected from radio to VHE 
gamma rays. In particular, it was detected by HESS, MAGIC and VERITAS showing a variable 
TeV emission on timescales of years although much faster variations, down to timescales of a 
day and less, have been observed [14, 15]. Cen A is at a distance of 3.8 Mpc. Its giant radio lobes, 
which subtend ~ 10° on the sky, are oriented primarily in the north-south direction. Also, Cen A 
has a jet with an axis subtending an angle to the line of sight estimated as 15° — 80° [0, ||]. 



1. Hadronic interactions 

We assume for protons a power law injection spectrum given by dNp/dEp = ApEp where 
ttp is the spectral index and Ap is the proportionality constant. 

We assume that the py interaction takes place when accelerated protons collide with target 



photons [|16p originated in the SSC process. Then, the energy loss rate due to pion production is 
given by t'p ^ = 1/(2 7^) j£ 5:t(£)<^ (e) £ /g~2yp'3fxx"^«(x) [|T^, where n{x) = dny/d£y{£y = x), 
<Jji{e) is the cross section for pion production for a photon with energy £ in the proton rest frame, 
1^ (£) is the average fraction of energy transfered to the pion, £0 = 0. 15 is the threshold energy and, 
7p = Ep/ni^p. The fraction of energy lost is /jjo py '^t'^/t'p y, (where t'^ ~ r^/Y is the expansion time 
scale) and the differential spectrum, dNy/dEy, of the photon-pions produced by p7 interaction is 
given by, 

"J 1 Tpohs ^ Tpobs 



dE J 
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where 



Ap,yoc5^"E^Ape-^yyil+z)-""E'^';'^ ^ L"^' df'" df^ (1.2) 
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and 



E"'^' , =5x10 



{\+z)H\-cosB) 



(1-3) 



where L is the observed luminosity, dt is the observed variability, Tyy is the optical depth y], 
5j) is the Doppler factor and is the distance to the source. 

On the other hand, for pp interactions we assume that the accelerated protons collide with 
thermal particles in the giants lobes. Then, the energy loss rate due to pion production is given by 

n'k,; 



t' 



'p '^PP^PP) 



[[Tq], where Opp = 30 mbarn is the nuclear interaction cross section, kp^ = 0.5 



pp 



is the inelasticity coefficient and n'^ is the comoving thermal particle density. The fraction of 
energy lost by pp is fjio^pp ~ t'd/^'pp ^^id differential spectrum of the photon-pions produced by 
pp interaction, dNy/dEy, is given by. 



dEJ 



obs 



pp-y 



(1.4) 



with. 



\2T^„ ^"P'^^ 



(1 



(1.5) 



where Hp is the thermal particle density, R is the distance to the lobes from the AGN core and T is 
the bulk Lorentz factor. The eqs. |0] and L4 represent the TeV energy range contribution to the 
spectrum. We have considered only one hadronic interaction at once. 

The input parameters used to fit the leptonic contribution to the SED of Cen A [§, |3|], M87 [§| 
and NGC1275[|8]] are shown in Table 1. In Table 2 the parameters used and their values to describe 
the TeV contribution to the spectrum considering py or pp interactions are shown. 



Parameters 


Cen A 


M87 


NGC1275 


A"*' (s) 


2.5 X 10* 


1.0 X 10^ 


3.17 X 10' 


L"*" (ergs-') 


5 X Ht^ 


1.0 X 10'*'' 


6 X lO'*^ 


6 (degrees) 


40 


10 


20 




1.47 


3.9 


2.9 


dz (pc) 


3.8 


16 


76 



Table 1. Input parameters used to fit the leptonic contribution of the SED 



Parameters 


Cen A 


M87 


NGC1275 


Parameters 


Cen A 


M87 


NGC1275 


dp 


2.82 


2.49 


2.7 


Up 


2.82 


2.49 


2.7 


E°*' (GeV) 


317.1 


674.9 


80,6 


np{cm'^) 


10-* 


10-1 


10-1 


Apj (MeV cm^ s)" 


1.37x10-- 


5.37 


4.37x10- 


App (MeV cm^ s)" 


5.9x10-" 


2.0 xlO-"* 


7.2 xlO-l"! 


Eo (TeV) 


1 


1 


1 


Eo (TeV) 


1 


1 


1 










R (kpc) 


100 


150 


200 



Table 2. Parameters used to describe the TeV spectrum considering p7 (left) or pp (right) interactions are shown. 
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2. Discussion and Conclusions 



When p7 interaction is considered and assuming that each photon-pion carries ~ 10%, we 
obtain a broken-power law (eq. 1.1) with break energies inside the energy range and spectral 
indexes observed by HESS, VERITAS and MAGIC detections of Cen A, M87 and NGC1275 
respectively. When pp interaction is considered and assuming that each pion carries ~ 18% and 
splits in two 7— rays, we also obtain energy ranges observed by the HESS, VERITAS and MAGIC. 
For Cen A the number density of thermal particles in the giants lobes and the distance lobes-core 
are known, Hardcastle et al. (2009). However for the M87 and NGC1275, these values are not well- 
determined. Fraija et al. (2012) shows the proportionality constant (App) as a function of distance 
to the core (R) for typical number density of thermal particles (Up), obtaining the minimum App 
values for this process. 

We observe that p7 interaction is much stronger that pp interaction. For Cen A, the extrapola- 
tion to ultra-high energies (UHE) of the proton spectrum could explain the observed UHE cosmic 
rays (UHECR) observed by the Auger Observatory only for case when pp interactions are the 
responsibles of the TeV emission . Otherwise the expected number of UHECR obtained is several 
orders of magnitude above the observations. We have taken into account the variability for M87, 
~ 10^ s, that corresponds to low state, for short timescale variations presented in flaring state this 
model cannot be accommodated. 
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